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Abstract 

We study a model of scalar field with a general non-minimal kinetic coupling 
to itself and to the curvature, and additional coupling to the Gauss Bonnet 
4-dimensional invariant. The model presents rich cosmological dynamics and 
some of its solutions are analyzed. A variety of scalar fields and potentials 
giving rise to power-law expansion have been found. The dynamical equation 
of state is studied for two cases, with and without free kinetic term . In both 
cases phenomenologically acceptable solutions have been found. Some solutions 
describe essentially dark energy behavior, and and some solutions contain the 
decelerated and accelerated phases. 

PACS 98.80.-k, 95.36+x, 04.50.kd 



1 Introduction 

The late time acceleration of the universe is one of the most important problems of 
modern cosmology, which is supported by astrophysical data from distant la super- 
novae observations [I], [2], [3], [I], cosmic microwave background anisotropy j5], and 
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large scale galaxy surveys [6]. The interpretation of astrophysical observations indi- 
cates that this accelerated expansion is due to some kind of negative-pressure form 
of matter known as dark energy ([7], [5], [H]). The combined analysis of cosmological 
observations also suggests that the universe is spatially flat, and consists of about 
~ 1/3 of dark matter, and ~ 2/3 of homogeneously distributed dark energy with neg- 
ative pressure. The dark energy may consist of cosmological constant, conventionally 
associated with the energy of the vacuum (TO], [TT], or alternatively, could came from 
a dynamical varying scalar field at late times which also account for the missing en- 
ergy density in the universe. A widely explored scalar field models are quintessence 
02], P2], HU, tachyon [IS], phantom [16], K-essence [T7J,[TH], and dilaton [19] (for a 
review see [7J). The scalar fields are allowed from several theories in particle physics 
and in multidimensional gravity, like Kaluza-Klein theory, String Theory or Super- 
gravity, in which the scalar field appears in a natural way. Another alternative to the 
explanation of the DE is represented by the scalar-tensor theories, which contain a 
direct coupling of the scalar field to the curvature, providing in principle a mechanism 
to evade the coincidence problem, and naturally allowing (in some cases) the crossing 
of the phantom barrier [20], [21]. From a pure geometrical point of view, the modified 
gravity theories, which are generalizations of the general relativity, have been widely 
considered to describe the early-time inflation and late-time acceleration, without the 
introduction of any other dark component, and represent an important alternative to 
explain the dark energy (for review see [22J and references therein). 
In the present work, we consider a model with non-minimal coupling to gravity 
[2"3l [2U [22], and with additional Gauss Bonnet coupling, specifically we focus in 
a scalar field model with kinetic term non-minimally coupled to gravity and to itself 
[26j [27J [28], with a new term containing the Gauss Bonnet (GB) 4-dimensional in- 
variant coupled to the scalar field, with the coupling given by an arbitrary function 
of the field. 

Despite the fact that the GB term is topologically invariant in four dimensions, 
and hence, by itself does not contribute to the equations of motion, nevertheless it 
affects cosmological dynamics when it is coupled to a dynamically evolving scalar 
field. Besides this, if the GB term is coupled to the scalar field through arbitrary 
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function /(</>), then this is the unique quadratic combination of the Riemann curva- 
ture tensor, that does not increase the differential order of the equations of motion 
(i.e. GB produces only terms which are second derivatives of the metric in the field 
equations). Therefore, the coupled GB term which preserves the theory ghost free, 
seems as a natural generalization of the scalar field with non-minimal kinetic coupling 
to curvature, which under certain condition is a second order scalar tensor theory. 
Both couplings my have origin in the low energy limit of higher dimensional theories. 
Thus, the kinetic coupling appears as part of the Weyl anomaly in N = 4 conformal 
supergravity J29J EO], while the GB coupling arises naturally in the leading order of 
the a' expansion of string theory (31], [32]. Besides that, the kinetic couplings to 
curvature are also present as quantum corrections to Brans-Dicke theory [33] and in 
different frames in modified gravity [22]. 

Some late time cosmological aspects of scalar field model with derivative couplings to 
curvature have been considered in [31], [35], [35]. On the other hand, the GB invariant 
coupled to scalar field have been extensively studied. In [37] the GB correction was 
proposed to study the dynamics of dark energy, where it was found that quintessence 
or phantom phase may occur in the late time universe. Accelerating cosmologies with 
GB correction in four and higher dimensions have been discussed in [33], |3j5], |40j . 
The modified GB theory applied to dark energy have been suggested in [41] . and 
different aspects of the modified GB model applied to late time acceleration, have 
been considered among others, in [22], [12], [13J, [11 . 

All these studies demonstrate that it is quite plausible that some scalar-tensor cou- 
plings predicted by the fundamental theory may become important at current, low- 
curvature universe. The aim of the present paper is the study of the possible effects 
of the Gauss-Bonnet coupling combined with the kinetic coupling, on dark energy 
cosmologies, in order to clarify whether they could be compatible with current phe- 
nomenology. 
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2 Field Equations 



Let us start with the action, for scalar field kinetic terms non-minimally coupled to 
curvature and coupled to Gauss Bonnet (GB) curvature 

(2.1) 



- V R pv (^(0)^0^0) - v(</>) + F 2 (<p)g 



+ s n 



where G is the 4- dimensional GB invariant Q = R 2 - AR^R^ + R^ pu R^ pa , S m is 
the dark matter action which describes a fluid with barotropic equation of state. The 
dimensionality of the coupling constants £ and rj depends on the type of function 
Fi(<f>), and the coupling F 2 (4>) is dimensionless. Besides the couplings of curvatures 
with kinetic terms, one may expect that the presence of GB coupling term may be 
relevant for the explanation of dark energy phenomena. 

Taking the variation of action (12. ip with respect to the metric, we obtain a general 
expression of the form 

- \g,uR = k 2 [T™ + T M J (2.2) 

where k 2 = 8nG, T™ v is the usual energy-momentum tensor for matter component, 
the tensor T^ v represents the variation of the terms which depend on the scalar field 
and can be written as 

T — qvt> _i_ t£ _i_ _i_ t gb (9 "i\ 

where T^ u , correspond to the variations of the standard minimally coupled terms, 
T^, TJj u come from the £ and rj kinetic couplings respectively, and Tff u B comes from 
the variation of the coupling with GB. Due to the kinetic coupling with curvature 
and the GB coupling, the quantities derived from this energy-momentum tensors will 
be considered as effective ones. The variations are given by 

It = V m0V,0 - ^V A 0V A - g, u V((j>) (2.4) 



1% =£ [ ( R ^ - \^R \ (F!(0)V A 0V A 0) + ^,V A V A (Fx(0)V 7 0V" 
- ^(V M V, + V„V„) (F!(0)V A 0V A 0) +R(F 1 (<p)V fl <pV u <P 



(2.5) 
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- \ (VaV^ (F!(0)V a 0V^) + V A V, (F 1 (0)V A 0V M 0)) 
+ ^V A V A (iM^V^V^) + ^V A V 7 (F 1 (0)V a 0V 7 



(2.6) 



and 

T^ B =4(jV p V,F 2 (0)]i? - v[VpVF 2 (^)]i? - 2[V"V M F 2 (0)]i2 I/p - 2[VV I ,F 2 (0)]f2 I/p 
+ 2[V p VF 2 {4>)]R pv + 2^[V' i V CT F 2 (0)] J R p(7 - 2[VV ff F 2 (0)]i2 /1/Wff ) 

(2.7) 

In this last expression the properties of the 4-dimensional GB invariant have been 
used (see [13], [37]). Variating with respect to the scalar field gives the equation of 
motion 



1 -.d, y=g {^RF x {cj>)dy + rjR^F 1 ( ( j ) )d u( p + &*</>)] + ^+ 



— - (tw^dy + V R^d^d^) - —tg = 



(2.f 



Considering the spatially-flat Friedmann- Robert son- Walker (FRW) metric, 

ds * = + a ( t )2 ( dr 2 + 



(2.9) 



And assuming an homogeneous time-depending scalar field <fi , the (00) and (11) 
components of the Eq. (|2.2|) . from 02.3H2.7p take the form (with the Hubble parameter 
H = a/a) 



H 2 = ^rPeff 



with p e ff given by 
1 



Peff 



^0 2 + V(<j>) + 9f# 2 Fi(0)0 2 + 3(2£ + //)///•,(•• 



-3(2£ + ri)HF 1 ( ( j ) ) ( j ) cf>--(2!;+ v )H 



dF 1 



-24# 



3^ ; 



and 



(2.10) 



(2.11) 



2# - 3F 2 = «V 



e// 



(2.12) 
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with p e ff given by 



Peff 



dp, . 

+ 4(£ + ^HF^c/xj) + 2(f + v)H-^<j) 3 




(2.13) 



where we have assumed scalar field dominance (i.e. = 0). The equation of motion 
for the scalar field (12.81) takes the form 



where the first three terms correspond to the minimally coupled field. In what follows 
we study the cosmological consequences of this equations, under some conditions that 
simplify the search for solutions. 

An important simplification of the Eqs. (I2.10H2.14)) takes place, under the restriction 
on £ and rj given by 



Under this restriction all the kinetic couplings that appear in the action (12 .ip become 
summarized in the term G^d^^d^cj), where G^ u = R pv — \g^ u R- In this case the field 
equations (l2.10H2.14"j) contain only second derivatives of the metric and the scalar field, 
avoiding problems with higher order derivatives [241 E3]- The modified Friedmann 
equations (I2.10p and (12. 12ft take the form 



...try ... . . 

+ 3H<P + — + 3(2£ + rj)HF x {4>)<j) + 3(14f + 5r))HHFi(<f>)<f> 




V + 2£ = 



(2.15) 




(2.16) 
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and 



-2H - 3H 2 =k 2 



lj>*-V(<f>)-([3H» + 2H)F 1 



2£if 2iq(0)# + 



+ 8# 



2 dF 2 



+ 8# 



+ IQHH—^S + 16i7 



The equation of motion reduces to 

4> + 3H<P+^ + 3£H 2 (^2F(0)0 + ^0 2 ) + 18£tf 3 F(0)0+ 



12£MF(0)0 - 24 ( M 2 + # 4 * 2 







(2.17) 



(2.18) 



Note that, independently of Fi ((/)), and for F 2 ((p) = const., assuming an asymptotic 
behavior of the scalar field as <fi = <po = const., gives rise to de Sitter solution, as 
can be seen from Eqs. (12TTU1) and (12TT41) . From ( 12~T0|) and ( l2TT4j) it follows that 



V = Vo = const and H = H = k^JVq/3. 

Next we try to study cosmological solutions of this model, giving rise to to accelerated 
expansion and acceptable behavior of the equation of state parameter (EoS). 



3 Power law solutions 

It is of interest to derive solutions that give rise to a power-law expansion (a ex t p ), 
as it's well known that this kind of expansion is characteristic of the evolution at the 
early stages of radiation and matter dominance, and also may respond for accelerated 
expansion. We will consider the effects of this new kinetic coupling in the cosmological 
dynamics, in the case of scalar field dominance (we further neglect any background 
radiation or matter contribution). 



3.1 Power law with free kinetic term 

In this section we will consider the model with the coupling Fi((p) = l/<f) 2 . This 
function has the important property of leaving dimensionless the coupling constant £ 
in fl2T6j) . 
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Replacing dF 2 /dt From (I2.16P into (I2.18p . and using -F\(0) = 1/0 2 , yields 



f |(0 2 ) + (5# 2 - + 3^ 3 | [ £ ) + 9^ + S^ 2 ^ 



+ #^ _ (> + ^ y + _3 # 2 1 /7 > „, = 



(3.1) 



where we have multiplied Eq. (1 2 . 1 8 [) by <fi and used (pd/d(j) = d/dt. Assuming the 
following time dependence for the scalar field and the Hubble function 

and integrating (13.1 p we find the potential 

9^-9^+5^-3 4 V(P-1) ^ ^ 1 {V - X)I2 C 

m = WTs) + WTYW^ + (3 ' 3) 

where we have replaced t = V2/4>. Solving the Eq. ( 12. 16ft with respect to F 2 , we find 

3& 2 + 3& + 1 / \ 1 1 2^+ 3 )/ 2 C 1 

2 1<W ~~ 4p 2 (p + 3) n V v^/ + l)/t 2 2 24p3(p + 3) 0p+ 3 1 ' J 

This solution gives rise to accelerated expansion, provided p > 1. A remarkable 
property of the potential (I3.3p . is that for the particular choice C = 0, becomes of 
the Higgs type. 

Let's consider the particular case of the model (12. ip without free kinetic term. In this 
case we are in the frames of the models ([15J H7J HH1 EH])- This also applies if the the 
slow-roll condition <fi 2 << V(<f>) is considered, which can take place in the contexts of 
dark energy or inflationary cosmology. 



3.2 Derivative coupling with V = 

Let us then, begin with the simple case of strictly non-minimal kinetic coupling, 
without potential term. In this simple case we may consider the general Eqs. (12.101) 
and (I2.14p . which after eliminating the free kinetic term an potential become 



H 2 = - 



9£H 2 F x (0)0 2 +3(2£+ V ) HF^-^+^Hj^F^^AH 3 ^-] (3.5) 



and 

3(2£ + v)HF 1 (<f>)4> 2 + 3(14f + hr])HHFi((j))(f) 2 + ~(2£ + rj)H±-{F x tf)+ 

o . , . jp (3-6) 

|(4£ + v)H 2 j t (F^) + 9(4£ + VW'F^ 2 - 24 (m 2 + if 4 ) ^ = 

where we have multiplied the Eq. ( I2.14p by <fi. Let's propose now a power law solution 
H = p/t, and the coupling F\ satisfying F\(f) 2 = C\. Then replacing in Eq. ( 13. 6 \ we 
obtain 

d Jl = Ct c = (4e + 277 - + 577)p + 3(4£ + y)p 2 ) C x 
dt 2 ' 2 8p(p - 1) 1 ' ' 

Integrating this equation gives 

F 2 (t) = l -C 2 t 2 (3.8) 

where we equated the integration constant to zero. Replacing this solution together 
with H = p/t and Fi<j) 2 = C% in Eq. (13. 5p . we obtain the following relation between 
C\ and p (setting k 2 = 0) 

' (3p- l)(2£ + r]-3£p-r}p 

giving rise to accelerated expansion provided p > 1. Thus for example taking F\ = 
F e 2 ^^ and limiting to those values of C\ > 0, the scalar field behaves as (from 
F^ 2 = d) 

= ^,]n[jSl) (3.10) 

and the GB coupling function F 2 from ( 13. 8 j) becomes 

F2 = 9^1 e ^ (3.11) 
2Gx 



(3p-l)£ 



Under the restriction (12.151) . the relation between Ci, C 2 and p becomes C\ 
and C 2 = 2 2 ~^ p 4 , which also gives rise to accelerated expansion for p > 1. For the 
same restriction ( 12.151) in the model without GB coupling, the power-law solution was 
limited to p = 2/3 
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3.3 Derivative coupling with potential 

From now on we will consider the restriction (|2.15|) . We will look for the shape of the 
potential V{<p) and couplings corresponding to the asymptotic behavior giving rise 
to power-law accelerated expansion. Dropping the free kinetic terms in ( 12 . 16[) and 
(I2.18P and replacing H = p/t we obtain 

4= v +4^-4^ <» 2 > 

and 

£ +4i M + -4 M - °4 M - - ° <»•»> 

where in this last equation we have multiplied by <fi and the first two terms of Eq. 
( I2.18P have been dropped due to the absence of the free kinetic term, using the addi- 
tional freedom provided by the GB coupling we can impose the following restriction 
on the potential 

dV _ Ap -Dgj, = (3.14) 
dt t 4 dt K J 

which converts Eq. (13.131) into a first order differential equation with respect to the 

variable {Fi<f) 2 ): 

4^M +i 4m- i 4m < 3 ' is > 

Integrating this equation with respect to (i<i0 2 ), gives 

t \ -(6p-4) 



Fi(j> 2 = Xo I - I (3.16) 

using the restriction ( 13. 14j) in Eq. (13. 12 j) . we get the following differential equation 
for the potential 



p - 1 dt kH 2 t 6 P 

Integrating this equation gives 



V = DT* + f -f'i + **.(?- Diftiri (3.18) 
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and from (I3.14p we find the expression for the GB coupling 



24p3(p + 3) (p+l)/t 2 (7p-3)(l-p) 

Finally, from Eq. (I3.16p . for a given i*\ we can find the scalar field as a function of 
time and then express the potential V and the coupling F2 in terms of the scalar field. 
But the expression (I3.16P also allows to propose the scalar field as function of time, 
and find the expression for the kinetic coupling Fx. Considering both possibilities, we 
first propose the kinetic coupling as Fx = F e~ 2 ^^ , which from (I3.16P gives (up to 
additive constant) 

4> = 3(p-l)(f> ln(P\ (3.20) 



\V 

making the rescaling 3(p — 1)0 O — > 0o an d replacing t in Eqs. (I3.18P and (I3.19P we 
find the following expressions for the potential and GB coupling 

V _ rf P~lJ P ~lWI>0 1 3(p - l)p 2 , 9 £Xo(P ~ 1)P 3 -2(3p-l)j>/0o /o 01 \ 

° + (P+I)« 2 tf + (3-7p)^ (3 - 21) 

and 

= Ctg +3 (p+3)^/^ _ V^o 2^/0 o , 3p 3 (3p - l)£xo*o -6(p-iMAto /o 22I 

2 24p3(p + 3) (p+1)k 2 (7p-3)(l-p) 1 ' J 

By assuming now for the scalar field the simple expression = t" 1 , from (I3.16P it 
follows g ^ 4 

Fx = = XotrY p - 8 (3.23) 

and the potential V and GB coupling F2 are obtained from (13. 18[) and (I3.19P by mak- 
ing the replacement t = 1/0. In this case the GB coupling and the scalar potential 
are combinations of different powers of the scalar field. This shows that in the case 
when the free kinetic term is neglected, we always can find the appropriate couplings 
and potential giving rise to power-law expansion. Note that looking at the time de- 
pendence in all considered solutions, the GB coupling decays faster than the coupling 
F 2 grows or decays (making the integration constant C = in (13. 4p and ( I3.19P ). 
Then the term F^Q decays with time, behaving properly in the current low-curvature 
universe. 
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4 Dynamically varying equation of state 



So far we have considered solutions of the power law type, which are characterized by 
a constant EoS parameter. Power law solutions may describe a decelerating universe 
with matter or radiation dominance, as well as an accelerating universe, but cannot 
describe the transition from one phase to another, due to the non-dynamical nature 
of the EoS parameter. In the present model one can explicitly construct a solution 
which admits the transition from the matter-dominant period to the acceleration 
phase dominated by dark energy. Considering as before the coupling Fi = 1/0 2 , the 
scalar field cf> = v2/t, and proposing the Hubble function 

ir(f) = 7 + f, (4.i) 

with 7 and p constants, the equation (13.1 j) takes the form (setting k 2 = 1) 
(pt 5 + 7 t 6 ) ^ + (p(l - p)t A - 2p 7 t 5 - 7 2 t 6 ) V + 3 7 4 t 6 + 12p 7 3 t 5 



(4.2) 



+ 37 2 (6p 2 -p + 3£ 7 2 ) t 4 + 6 7 (2p 3 - p 2 - £ 7 2 + 6p^ 2 ) t 3 
+ (3p 4 - 3p 3 + 5 7 2 - 21p£ 7 2 + 54p 2 £ 7 2 ) t 2 + 2 7 (l8p 3 £ -hp-2)t 
+ 9£p 4 - %p 3 + 5p 2 - 3p = 
To solve this equation we propose a particular solution of the form 

nt) = V + j + § + % + % (4.3) 
Replacing this solution in ( 14. 2 p and solving, we find two solutions: 

V a = , A = 0, B=^~ l \ £7 = 0, D= *& + 7 = 

p+ 1 p + 3 

(4-4) 

which reproduces the previous result for power law expansion (see ( 13 .3p with C = 0), 
and 

T/ Q 2 c QP + 5P 2 +P 3 ~1 2 . _ 

^0 = 37, e = — 2(p _ 6h2 > ^ = 6 ^ 

B _ 3 (3 7 2 + 6p - Zip 2 - p 3 ) c _ 3 (3p 4 + lip 3 - 6p 2 - 3p 7 2 + 4 7 2 ) g 
2(p — 6) (p — 6)7 

_ 12 7 2 - 26p7 2 + 9p 2 7 2 + 18p 3 - 9p 4 - 9p 5 
~ 2(p - 6)7 2 
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w = -l + (U 



The potential can be written explicitly in terms of the scalar field, using t = y/2/<j> 
V(<P) = Vo + fy+ f 2 + -^f + j0 4 (4.6) 

replacing V from (14.31) in Eq. (12.161) . using (14. 5 p and solving for the GB coupling F 2 , 
we find 

2p(6-p) 7 2 (p0 + v / 27) (4 7) 

(3p 4 + 9p 3 + Qp 2 + (4 - 3p)7 2 ) (In a/2 - \n( P (f) + y/2'y)) 
8p 2 (6 — p)7 2 

where after solving with respect to time variable, we have replaced t by a/2/0. This 
expressions are valid for any p ^ 6, and therefore the model may reproduce known 
cosmological scenarios. As follows from (14. ip the EoS parameter is 

2 p 
3~7t 2 + 2-ypt + p 2 

for small t we can approximate w ~ — 1 + 2/(3p), and the second term in (14. ip dom- 
inates, giving a scale factor a ~ t p . Therefore, if p = 2/3 or j) = 1/2 the matter 
dominated or radiation dominated period could be realized. At late times, the first 
term 7 in (14. ip dominates and the universe evolves towards asymptotically de Sitter 
space time. Hence, the model is consistent with current phenomenology. 



Let's consider now the model without free kinetic term, and kinetic coupling 
Fi{4>) = 1/0 2 . Working with the variable x = In a, the Eqs. (I2.16P and (I2.18P (after 
multiplying by </>) take the form 

** = V + B eHV-KH*™ X =\f (4.9) 
k z dx <p dx 

and 

d JL + 9£H 2 ^ X 2 + 3£# 4 ^ + 18£#V - 12H 2 (h + H 2 )^ = (4.10) 

dx dx dx \ J dx 

where we used 4-.—R4-- 

at ax 

Considering the particular form of the scalar field = <poe~ ax and the GB coupling 
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F% = Aln(0/0o); X takes the value x = ~ a an d dF2/dx = \\ = — Aa. Then, taking 
the derivative of Eq. ( 14. 9 p and replacing dV/dx in ( I4.10p we obtain the following 
evolution equation for the Hubble parameter 



1 "' hh - 3a (£a + 4A) H 2 ^ + 2a (3£a + 4A) H A 



k 2 dx ^ x 
the solution to this equation is given by (setting k 2 

1 







(4.11) 



1) 



H\x) = - 3a(£a + 4\)W 



2(3ag+4A) , 



_g 3(a£+4A) 3a(£a+4A) 



(4.12) 



3a« + 4A) 

where W is the Lambert W function and C is the integration constant. From (I4.12p 
and the properties of the W function, it follows that for a and £ satisfying the 
inequality a(a£ + 4A) < 0, H 2 will always positive (otherwise we may enter in a 
region where W becomes complex, see [26] for details). Note that for A = we 
obtain the solution described in [26], and A = — 3ct£/4 gives rise to de Sitter solution. 
Choosing A = 1 and in the redshift variable z, the Hubble parameter takes the form 

/ r c -i \ -i 

H 2 (x) 



g 3a(£a+4) 2(3ag+4) 
'1 + Z \ 3(aS+4) 



3a« + 4) 



(4.13) 



and evaluating the equation of state (EoS) parameter w e ff(z) 
(14. 13p one obtains 

8 + 6r] 



1 + from 



-1 - 



9K + 4) 



1 + W 



e 3a(^+4) , 

3a(?7+4) v x T ^/ 



2(3.;+4) 



(4.14) 



where 77 = a£. In order to keep i7 2 positive, a and 77 should satisfy a(rj + 4) < 0. 
From the expression (I4.14p it follows that the behavior of w at the limits z — > 00 
and z — > — 1 does not depend on the values of a and C. Thus, a and C can be 
used to have an appropriate value of the current EoS parameter wq. An interesting 
value for 77 is 77 = —28/9, which gives a radiation dominated universe at very early 
times w\ =1/3 and cosmo logical constant solution at future to , = — 1, 
independently of the values of a and C . Taking a = — 1 and C = 200 the current 
value of the EoS is wq ~ —0.98. Another behavior takes place for 77 = —24/9. In this 
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case the EoS describes an effective quintessence evolution, as it starts in the lower limit 
of the accelerated phase w\ = — 1/3 and ends in a de Sitter phase w\ = —1. 
And the current value of the EoS is Wq ~ —0.97 with a = — 1 and C = 100. An 
universe evolving from a matter dominated phase at high redshifts w\ = 0, is 
obtained for rj = —44/15. This solution evolves towards de Sitter phase at far future, 
with an EoS at the present Wo ~ —0.99 (for a = —0.1 and C = 50). There is also a 
possibility for a universe in the phantom phase, with an EoS evolving from w = —1 
at 2 - > oo to w < - 1 at z - y — 1 (i.e. starting in a de Sitter phase and ending in 
phantom phase) with an EoS currently less than —1. Thus, for instance r\ = 1, gives 
w = —1 at z — y oo and w = —59/45 at z — y —1, with the current value Wq = —31/24 
for a = —1 and C = 1, or taking C = 50 we have Wq ~ —1.17. In this case the model 
describes pure dark energy behavior. 

From the expression for H 2 ( 14.13)) . it follows that the effective energy density given by 
p e ff = 3H 2 / k 2 shows two types of behaviors: for a solution starting from radiation or 
matter dominance, or purely quintessence behavior (i.e. rj = —28/9,-44/15,-24/9), 
the evolution begins from Big Bang singularity, as p e ff\ = oo and presents no 
future singularities, p e ff \ = . And for solutions describing a phantom phase 
(i.e. rj = 1), p e ff tends to zero at high redshifts, validating the assumption that 
the dark energy was negligible at early times, but presents asymptotical singularity 
(p e ff\ = oo), characteristic of the phantom behavior. 

Another interesting solution takes place if we impose the restriction £ct + 4A = 0, in 
Eq. ( 14. lip , which after integration gives the result (setting k 2 = 1) 

H 2 = = (4.15) 

4£a 2 x + C 4£a 2 lna + C V ; 

where C is the integration constant and we replaced x = In a. Integrating this 
equation with respect to time, gives the scale factor 

1/3 



a(t) = ao exp 
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f 2/3 



(4.16) 



hence, when t is small we may expand (14.161) keeping the first two terms, and the 
scale factor behaves as a ~ t 2//3 . Therefore, the matter-dominated period could be 
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realized. On the other hand, when t is large, from (14 . 1 6 j) the Hubble function becomes: 
H = ^7173 1 ■ Then w e ff = — 1 + 2 ^ 6 ^ 2 /l ' ? an d therefore the universe accelerates 
and is asymptotically de Sitter space. 

5 Discussion 

We have considered a model of scalar field that contains non-minimally kinetic terms 
coupled to curvature and the Gauss Bonnet four dimensional invariant coupled to the 
scalar field. This GB coupling has the advantage that does not make contributions 
higher than second order (in the metric) to the equations of motion, and therefore 
does not introduce ghost terms into the theory. Nevertheless, this new coupling in- 
troduces a new degree of freedom that could be exploited to enrich the number of 
possible evolutionary scenarios that explain the nature of the dark energy. Power 
law solutions have been found for different configurations of the model, namely, the 
complete model, and the model without free kinetic term, with and without potential. 
Under the restriction 77+ 2£ = 0, and neglecting the free kinetic term, the GB coupling 
allows power-law solution with accelerated expansion (p > 1), improving the results 
in absence of GB coupling where the only decelerated expansion with fixed p = 2/3 
was possible [26]. The power-law solutions with (p > 1) are also possible in absence 
of potential, or with potential and GB coupling as a combination of exponentials of 
the scalar field. 

Studying cosmo logical solutions with varying EoS, we considered two cases. In the 
complete model we found an interesting solution that contains the two asymptotic 
behaviors corresponding to early time (matter or radiation dominated universe) and 
late time dark energy dominated universe, evolving towards the de Sitter phase |22j . 
In absence of free kinetic term, we considered the couplings F\ = 1/0 2 (which leaves 
the coupling constant £ dimensionless) , F2 = A ln(0/0o)> and assumed the scalar field 
= 0o e_QX (which is a reasonable choice as e~ ax = (1 + z) a ), we obtained a vari- 
ety of solutions with different asymptotical desirable behaviors, and that reproduce 
the current observed value of w. The solutions with rj = —44/15,-28/9 describe 
an universe starting asymptotically in matter or radiation dominance, and evolving 
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towards the future de Sitter solution, with current values wq ~ —0.99, —0.98 respec- 
tively. Therefore this solutions present transition between decelerated and accelerated 
phases. Taking for instance 77 = —24/9 gives rise to an universe in quintessential phase 
starting with w = —1/3, with current Wo ~ —0.97 and ending in a de Sitter phase 
with w = — 1. There is also a possibility for a solution with universe in essentially 
phantom phase for 77 = 1 beginning in de Sitter phase, with current Wq = —31/24, 
and evolving towards a state with w = —59/45 (note that in all cases the current 
Wo may vary, depending on the choice of a and C). This last two types of solutions 
belong to purely dark energy description and need the addition of matter content 
to achieve a phase transition between decelerated and accelerated phases. All these 
results are based on the properties of the Lambert W function, which varies very 
slowly in an infinite (redshift) interval. These results, that were obtained thanks to 
the GB coupling, generalizes the one obtained in [26] as in that case the argument 
of the Lambert W function depended on fixed power of the redshift ((1 + z) 2 ) giving 
only one possible cosmological scenario. 

A remarkable property of the solution (14.121) is that the cosmic evolution was obtained 
(i.e. was not proposed a priori) as a result of solving the dynamical equations for given 
couplings Fi((p), F 2 ((p) and scalar field (instead of the potential V), reflecting in 
some sense predictive power of the model, as the obtained solution is compatible with 
current DE phenomenology. In the DE models, mainly due to the complexity of the 
involved equations, one usually gives a dynamical evolution a priori, and reconstruct 
the model satisfying this solution. Is important to note, that in all solutions con- 
sidered here, the Gauss-Bonnet invariant Q decays faster than the coupling F 2 (<f)) 
grows or decays, so that at late times the term F2{<p)Q decays and therefore, does not 
substantially affect the current restrictions on the time variation of the gravitational 
coupling with respect to that evaluated with the kinetic coupling 12"T] . 
In conclusion, we have shown that the combined effect of non minimal kinetic cou- 
pling to curvature and the Gauss Bonnet coupling to scalar field, could represent an 
interesting source of dark energy, which might play an important role in the explana- 
tion of the current (late time) cosmological dynamics. We have concentrated mainly 
on the scalar kinetic coupling coupling Fi = 1/0 2 , while the study of other type of 
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couplings and their role in the late time cosmology will be considered elsewhere. 
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